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Abstract. A few examples are given of Ctiandra's work on statistical and stochastic problems 
that relate to open questions in astrophysics, in particular his theory of dynamical relaxation in 
systems with inverse-square interparticle forces. The roles of chaos and integrability in this theory 
require clarification, especially for systems having a dominant central mass. After this prelude, 
a hypothetical form of bosonic dark matter with a simple but nontrivial statistical mechanics is 
discussed. This makes for a number of eminently falsifiable predictions, including some exotic 
consequences for dynamical friction. 
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1. Introduction 

Probabilistic and statistical considerations dominated several of Chandrasekhar's contri- 
butions to astronomy and astrophysics. Among his less-well-remembered contributions 
to astronomy proper, as opposed to astrophysics, are a series of papers with Miinch on 
brightness fluctuations in the star fields of the Milky Way [beginning withT] and one on 
the inference of the distribution of stellar spins from observed values of v s'mi \2\. In a 
paper that deserves to be better remembered lH], Chandra introduced the use of structure 
functions to describe the refractive-index fluctuations that degrade ground-based astro- 
nomical images ("seeing"). These methods have been used ever since by students of 
atmospheric and interstellar scintillation, although usually without recognizing the debt 
owed to Chandrasekhar. Better known is Chandra's magisterial review article on random 
walks and other stochastic processes [41 . 

The problems of this type for which Chandra is most famous are those of dynamical re- 
laxation and dynamical friction. He was not solely responsible for these developments — 
the problem goes back at least as far as Jeans [5J, and a fairly complete solution for the 
velocity diffusion coefficients of a plasma subject to Coulomb encounters was given in 
the 1930s by Landau [61, but Chandra's subsequent emphasis on the probabilistic aspects 
was distinctive. He considered the probability distribution of the gravitational force expe- 
rienced at a given point in space due to an infinite field of point masses whose locations 
form a Poisson process (the Holtzmark distribution), for which Prob(> /) oc /^^^^ at 
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large values of the magnitude / of the force, and he then considered the effect of the mo- 
tion of these gravitating masses on a test mass |7 1. Rather than describe these as a series of 
impulsive two-body encounters, as done in most textbooks [e.g. 8 1 and indeed in his own 
monograph ||9l, Chandra treated the summed gravitational forces of the field particles as 
a time-dependent stochastic process. In Chandra's hands, this approach described the dif- 
fusion of the test particle's velocity, but not (at least in the first instance) the back-reaction 
on the field particles; therefore it predicted indefinite increase in the test particle's kinetic 
energy. Chandra recognized the need for a drag term — dynamical friction — to preserve a 
Maxwellian equilibrium distribution ifTOl . 

It is well known that the usual treatments encounter a logarithmic divergence, the 
usual cure for which is to estimate lower and upper bounds for the impact parame- 
ters of two-body encounters. The divergence at small impact parameters is straightfor- 
ward to remove by replacing the impulse approximation with exact formulae for hyper- 
bolic close encounters, as Chandra demonstrated, so that it suffices in most cases to set 
&min ^ G'(mt + mf these symbols having the usual meanings. But the divergence 

at large impact parameters is more problematic. It is conventional to remark that the treat- 
ment of the field particles as statistically uniform breaks down at distances comparable to 
the size or density scale height of a self-gravitating system (R) and therefore to cut off 
the impact parameter at 6,„ax ^ R, so that (since {v^^i) ~ GnifN/R) the argument of 
the logarithm is of order N, the total number of particles. Chandra, however, argued for a 
cutoff at the mean interparticle distance, D ~ R/N^/"^. This choice seems to be in some 
tension with his view of the gravitational field as a collective process. 

Although it is now accepted that fomax/&min oc in self-gravitating systems, much 
ink has been shed over the proper choice of the constant of proportionality. As a prac- 
tical matter, it makes little difference for large A^. Moreover, having worked in stellar 
dynamics for his Ph.D. and several years thereafter, the speaker and first author of this 
talk (JG) takes the view that the whole discussion is futile because there can be no such 
thing as a universally correct choice. The logarithm must depend at some level upon the 
details of the density profile of the system and upon where the test particle orbits within it. 
Furthermore, the root of the difficulty with large impact parameters lies in the insistence 
upon Markovian approximation; that is, the assumption that the probability distributions 
or moments of the test particle's velocity and position can be described by differential 
equations. This may be a very good approximation for gas particles interacting via short- 
range collisions, or even for Coulomb encounters in a plasma, since they are naturally 
cut off at the Debye length. But gravitational encounters have a duration that increases, 
on average, in proportion to impact parameter, b. Even in an infinitely system of field 
particles, the evolution of the probability distribution P[v, t\vf), io] for the test particle's 
velocity, should be insensitive to encounters more distant than b{t, to) = \t — ioK^^rei)^^^' 
stars more distant than this exert a net acceleration that is effectively constant during the 
interval t ~ to. The divergence of the net acceleration as the size of the system tends to 
infinity is removable (by shifting to a freely-falling frame) because it is constant for all 
finite time intervals. To reformulate dynamical relaxation along these lines would require 
integral rather than differential equations in time, however The slight theoretical ambi- 
guity represented by the usual logarithm is a price well worth paying for convenience in 
practical applications. 

There is another bit of unfinished business regarding dynamical friction that is proba- 
bly more urgent, especially because of the discovery of extrasolar planetary systems and 
the consequent revival of interest in celestial mechanics. (Chandra did not work in ce- 
lestial mechanics. Perhaps it seemed to him to be a closed subject. But in view of the 
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mathematical elegance of the subject and the new questions brought to light by computer 
calculations as well as planetary discoveries, it is likely that he would take interest if he 
were alive today.) This unfinished business concerns the relationship between dynami- 
cal relaxation and chaos. Chandra's analysis presumes that two given particles encounter 
one another just once since he treats their unperturbed trajectories as rectilinear This 
eliminates the possibility of resonances, since resonances require repeated encounters. 
Celestial mechanics, however, is rife with resonance; and chaos begins where resonances 
overlap 1 1 1 1. It seems that one uses the language of celestial mechanics when the number 
of gravitating bodies is small, N < 10, but that of dynamical relaxation when N > 10'^. 
Ambiguous intermediate cases exist: For example, the minimum number of particles in 
the core of a cluster of equal point masses undergoing core collapse (and rebound) is 
thought to be < 30 1 12 1; and the number of significant bodies in a young planetary sys- 
tem might well lie in the range 10-100 before this number is reduced by consolidation or 
ejection. It would be interesting, and perhaps useful for application to such intermediate 
cases, to understand these two domains of relaxation and resonant interaction as different 
limits of a unified statistical theory of Nbody systems. Some steps in this direction have 
been taken lfT3HT4ll . 

After this prelude, the remainder of this talks is devoted to a rather different topic: a 
(probably counterfactual) species of interacting dark matter. The connection to Chandra's 
work is admittedly tenuous, and consists mainly in the unusual properties of dynamical 
relaxation for such a dark matter. 

2. Dark matter and its possible collisionality 

Nonbaryonic dark matter (hereafter DM) makes up some 22% of the critical density 
SHq/SttG ifTSl . Notwithstanding experimental searches and suggested evidence for anni- 
hilation signals in cosmic backgrounds, DM is so far known only through its gravitational 
effects, including galactic rotation curves, virial temperatures of X-ray emitting gas in 
clusters of galaxies, and the power spectrum of temperature fluctuations of the cosmic 
microwave background. In compliance with Occam's Razor, DM is therefore usually 
presumed to have a minimum of non-gravitational interactions. However, were he alive 
today, William of Occam would presumably allow DM to have a minimum of additional 
properties if those would resolve important problems in particle physics as well as astro- 
physics and cosmology. In this spirit, the axion and supersymmetric weakly interacting 
massive particles (WIMPs) have been put forward as dark-matter candidates. The ex- 
istence of axions might explain CP conservation in the strong interactions |16|, while 
supersymmetry would unify fermions and bosons, and perhaps also the strong, weak, and 
electromagnetic coupling constants at high energies ITtI . WIMPs are favored by the con- 
sideration that if they were produced thermally in the early universe, then the present-day 
density of DM is consistent with WIMP masses in the range 1-100 GeV and annihilation 
cross-sections comparable to what would be expected from the weak interactions ifTSl . 

Effectively non-interacting DM would be consistent with most of what is known or 
believed about the present universe, and so it is conventional to treat DM as collisionless. 
This assumption lends itself to predicting the dynamics of DM via Nbody simulations. 
Occasionally, however, attention is drawn to apparent discrepancies between collision- 
less, dynamically "cold" DM and certain features of galactic structure on small scales. 
It is suggested on this basis that the DM may have additional properties such as an ap- 
preciable initial velocity dispersion ("warm" DM), or even some collisionality. One such 
discrepancy is the dearth of small satellite galaxies compared to the abundance of DM 
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subhalos predicted to lie within the DM halos of luminous (> L*) galaxies like the Milky 
Way lfT9l . This, however, might be explained by a low efficiency of star formation in 
the subhaloes |20|. Another discrepancy is the prevalence and persistence of apparently 
rapidly rotating galactic bars despite the dynamical friction exerted on these bars by the 
DM halo|21 1. But perhaps the dynamical friction is reduced by resonant effects too del- 
icate for standard Nbody simulation methods ll22l . A decade ago, the most talked-about 
discrepancy was an apparent minimum core radius for dark halos inferred from rotation 
curves of dwarf galaxies Il23lj26l . To solve these apparent problems, Spergel & Steinhardt 
||27| suggested that the (WIMPy) DM might have a coUisional mean-free-path in the range 
1 kpc < Amfp < 1 Mpc. This would efficiently erode DM subhalos and cause cuspy DM 
cores with profiles p^M c< r~" in the range < a < 2 to be smoothed out, at least 
initially, because the dynamical temperature associated with such profiles increases with 
radius. However, as the authors acknowledged, the same process of thermal conduction 
in a self-gravitating system leads eventually to core collapse, i.e. an indefinite increase in 
central density. They tried to argue that the timescale for collapse exceeds a Hubble time. 
But the timescale for core collapse cannot greatly exceed that for cusp smoothing except 
in carefully engineered profiles that would not prevail generally. 

Galaxy clusters in collision provide perhaps the most decisive limits on self-interacting 
DM. Comparison of the distributions of gas (traced by X-rays) and dark matter (traced by 
weak gravitational lensing) in IE 0657-56, the "Bullet Cluster," has been used to place 
a limit on the DM scattering cross section per unit mass: a/m < 1.25 cm'^ g^^ at 68% 
confidence 1281 . This translates to Amfp > 0.5 Mpc for the DM in the Solar Neighbor- 
hood. 

3. Repulsive dark matter: Overview 

A different sort of self-interacting DM was proposed at about the same time as ll27l : 
bosonic, like the axion, but with a repulsive interaction Il29ll30l . Also like the axion, these 
particles were supposed to be born in a Bose-Einstein condensate in the early universe. 
The interaction would naturally be short-range, in fact pointlike, if it corresponded to 
the non-relativistic limit of a massive scalar field with a momentum-independent self- 
interaction term V{(l)) > 0. In the condensate, pressure would vary only with density. 
In fact if V{(j)) = K(j>'^, with k a dimensionless coupling, then it was argued by fSO) 
that in the non-relativistic limit P <C pc^, the pressure and density would be related as 
in an n = 1 Emden polytrope, P — Kp^, with K — 3K/2m'^ in units fi = c = 1, 
where m is the particle mass. As Chandra himself discussed in his monograph on stellar 
structure Oli . the cases n G {0, 1, 5} are the only polytropes for which a (nonsingular) 
self-gravitating spherical equihbrium can be found analytically; for n = 1, the density 
profile is p{r)/p{0) = sm{r / a) / {r / a), < r < 7ra, with scale length a = ^jKj^-nG. 
If all of the DM were still in the condensate, then all nonrotating dark halos in virial 
equilibrium would have this size and density profile independent of their masses. Clearly 
this would not be acceptable, especially if the scale length a ^ 1 kpc as required to match 
the DM cores of dwarf galaxies. More complex and extended profiles can be obtained 
if the DM has a finite temperature, as discussed below, or where it has not yet reached 
thermal or even virial equilibrium. 

It was argued by t29 J and |.30i that if present-day repulsive dark matter (hereafter RDM) 
derives from a relativistic scalar field, then it has acceptable behavior in the early universe, 
i.e. it does not suppress large-scale structure or unduly affect primordial nucleosynthesis. 
There is some doubt whether even a massive scalar field has a sensible non-relativistic 
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limit, but we shall not worry about that here. We treat the RDM as an assemblage of 
nonrelativistic point particles of mass m having a two-body interaction potential whose 
range is small compared to the particles' Compton wavelength, so that U{xi — X2) — > 
Uo6{xi — X2)- The constant Uq = U{0) is the fourier transform of U{Ax) at zero 
momentum. If the relativistic correspondence were valid and V{(l)) = K(j>'^, then Uq = 
SnfV^ /rn?c at tree level. We further assume that that the particles are conserved. In the 
relativistic correspondence, this would require either that k be very small, or else that the 
particles be protected from annihilation by some symmetry, e.g. by making (j) a charged 
scalar field and its action invariant to global phase changes, since otherwise they would 
have an annihilation cross section at 0{k^ /m?) (annihilation of two particles into one is 
kinematically forbidden, so the lowest-order graph has two interaction vertices). 

In terms of the two basic parameters of the nonrelativistic model, the core scale length 
is a = y^C/oTi^rGrrz?, whereas the scattering cross section per unit mass is, to lowest 
order, a/m = itiUq /nh'^. Since these involve independent combinations of m and Uq, 
the degree of collisionality of RDM is independent of its minimum core size. So it seems 
possible to evade the constraint on CTscatt/w set by the Bullet Cluster [28 1. Also, since 
Cann/cscatt oc /t^ oc (m'^Uo)'^, reducing the collisionality reduces the annihilation rate 
faster. As will be shown, however, a long mean-free path leads in this species of DM to 
unacceptable halo density profiles, at least if collisions are frequent enough for thermal 
equilibration over a Hubble time. 

4. Superfluidity, dynamical friction, and galactic bars 

The repulsion makes the condensate, a large number of bosons in the same single-particle 
state, a superfluid. Superfluidity involves a penalty for excitations out of the degenerate 
state, which arises from the repulsive interactions. As a consequence of their indistin- 
guishability, the repulsion between a pair of RDM bosons in distinct momentum states 
\p) & \p') [{p\p') = 0] is greater than it would be if both were in |p) by [e.g.l32l 

AE = iyU{p-p') + ~ ^' ' p. + ^\^'\ (1) 

since the exchange energy reduces to i'U{0) for a short-range potential. Consequently, a 
rigid obstacle ("spoon") passing through the condensate at relative velocity v causes no 
dissipation if It^l < Wcrit = \/2i\Jrn because, viewing the problem in the rest frame of 
the spoon, the energy liberated by reducing the kinetic energy of a particle to zero is less 
than the exchange energy. Similarly, two streams of RDM condensate pass freely through 
each other if their relative velocity is less than v^t'crit even when the mean-free-path 
would otherwise be short. Here i/^ is the total number density of the condensate summed 
over both streams. 

At the time of |30|, one of us (JG) thought that the superfluidity of RDM implies 
that there should be no dynamical on a moving potential, such as a galactic bar, insofar 
as the core of the halo is dominated by the condensate and Wbai < ^^ciit- Our current 
view is somewhat different. In the limit that the mean-free-path is large compared to 
the dimensions of the bar or the core, the original argument is perhaps correct. In the 
opposite limit, while it remains true that the dynamical friction vanishes below a critical 
speed ^ Wcrit, the reason is not because the RDM is a superfluid, but rather because it is 
an ideal fluid in the classical sense, i.e. one in which the viscosity based on the mean- 
free-path and thermal velocity is negligible. The important difference between a galactic 
bar interacting with DM gravitationally and a spoon interacting with a laboratory fluid 
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is that the gravitational potential is smooth, whereas the spoon has a surface at which 
a no-slip boundary condition applies to the normal (non-superfluid) component. Thus, 
even in laminar flow, the spoon exerts a viscous force on the normal component that 
vanishes more slowly than linearly with the viscosity, probably as Re^^l"^, because there 
is a laminar boundary layer. At moderately high Reynolds number, i?e, the boundary layer 
on the rear side of the spoon becomes unstable and triggers turbulence. For a smooth 
large-scale potential like that of a rotating bar, however, there is no such surface and 
no such boundary layer, and therefore perhaps no turbulence. This absence of drag in 
the limit of small Amtp would seem to hold even if the RDM had a substantial normal 
component, except insofar as turbulence may arise even without boundary layers in a 
smooth velocity field at large Re. 

On the other hand, the bar or other moving potential may experience a wave drag if it 
couples to a waves whose phase velocity, whether linear or angular, matches that of the 
potential itself. In a pure condensate, below Udit, the only significant waves are sound 
waves at speed Cg = \/2Kp = Vcnt/V^, so that a linearly moving potential experiences 
no wave drag if its velocity is subsonic, though a rotating bar could excite sound waves 
at large radii were \r X ^hax\ > Cs- One can see this explicitly from the quantum- 
mechanical equation of motion for the condensate wave function (Gross-Pitaevskii or 
nonlinear Schrodinger Eqn.): 

ih— = - — V^* + C/o|*P* + ($.cif + *cxt)* , (2) 
ot 2m 

in which <i>oxt & 'i'scif are the gravitational potentials of the bar or other external per- 
turber and of the condensate itself; $scif satisfies Poisson's equation with the mass den- 
sity pc = mj'I'p of the condensate as source. Adopting the usual Jeans swindle, consider 
perturbations to a background state of uniform density and vanishing <!>. The unperturbed 
wavefunction has a constant modulus = but a time-dependent phase, be- 
cause the term involving C/q in eq- ^ doesn't vanish in the background state. One sets 
'J' — ^o(t)[l + e(r, t)] and expands eq. (|2]l to first order in the real and imaginary parts 
of e, treating $cxt and $soif as well as e as first-order quantities. In the unforced case 
^ext = 0, the dispersion relation for Fourier modes £{r,t) oc exp{ik ■ r — iojt) becomes 



w| = (2C/ofc2 - 47rGm)i.o + . (3) 

As usual with the GP Equation, the last term on the right represents single-particle exi- 
tations; it is small for long-wavelength modes k <C y/2mcs/h, with sound speed Cg = 
^/2UqI'o. If one neglects this term, then eq. (|3]l matches the results of a Jeans analysis for 
a classical ideal fluid, oj^ = Cg fc^ — AnGp^) [e.g. 8 1. If one now introduces a rigid perturb- 
ing potential that moves at constant velocity V through the condensate, <I'oxt(''" — Vt), 
and writes "I>cxt (fc) for the spatial fourier transform of this potential at any time, then after 
transients have decayed, the component of the wave drag along V is 

i^drag --yjj^ Si^'' k . y)|fe$ext(fc)p. (4) 

Since the fourier components ^cxtC^) are negligible at fc > y/2mcs/h, this is effectively 
the same drag as for an ideal fluid with equation of state P = Kp^. The quantum- 
mechanical nature of the condensate plays no direct role. To the extent that the self-gravity 
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of the RDM is slight on the scale of the perturber, in other words $oxt(fc) is unimportant 
for < AnGpo/c^, the drag vanishes unless V is supersonic. 

The formal result (|4]i for the wave drag clearly holds more generally for ideal fluids 
in which the dispersion relation may differ from ([3]). Thus it should hold even when the 
RDM has a normal (nondegenerate) component, as it must at finite temperature. In a 
realistic case where the background state is not uniform, however, the RDM gas would be 
stratified, i.e. it would have entropy gradient parallel to the background gravitational field, 
so that waves restored by buoyancy (internal waves/ g modes) might be excited at subsonic 
velocities. As will be shown in f|6] an isothermal RDM halo at nonzero temperature 
consists of a core that is almost pure condensate, and has nearly the n — 1 Emden profile, 
surrounded by an extended nondegenerate "atmosphere," with a sharp shelf in the density 
profile at the edge of the core. Although we have not calculated it, the coupling of a 
subsonicly rotating bar potential to the g modes would probably occur mainly at this shelf, 
with a drag proportional to the density in the nondegenerate component there; because 
that density is much less than the central density of the degenerate core, the drag on the 
bar would be expected to be considerably reduced compared to the estimates of [21] for 
collisionless DM at the same overall mass-to-light ratio. 

5. Thermodynamics 

Apart from the above clarification of its implications for dynamical friction, our main 
progress concerning RDM since [30 1 has been to construct isothermal spherical halo mod- 
els, which has lead us to a lower bound on the colUsionality of the RDM that is, while 
somwhat imprecise, nevertheless orders of magnitude less than what has been inferred 
from the Bullet Cluster — hence ruling out RDM, we think. 

For this purpose we needed to derive the equation of state of RDM at finite temperature, 
which we did by calculating the partition function of RDM considered as an almost ideal 
gas. "Almost" means that the scattering length Agcatt = mUo/h^ ^TrcTgcatt) is much 
less than the de Broglie wavelength AdB = h/\/2mkT. Then, just as in a classical 
dilute gas, the (non-condensate) particles can be regarded as belonging to well-defined 
momentum states between collisions. On the other hand, AdB should be small compared 
to galactic scales so that we may consider the partition function of a homogenous boxful 
of RDM. 

The almost-ideal boson gas is a well-studied physical model for a superfluid. The 
standard treatment in textbooks [e.g.[33l emphasizes the dynamics of quasi-particle exci- 
tations when most of the particles are in the condensate, and is couched in terms of anni- 
hilation and creation operators for these excitations. Our approach is cruder but, we hope, 
sufficient to determine the pressure in local thermodynamic equlibrium (LTE). Insofar as 
possible, we follow elementary treatments of a noninteracting ideal boson gas, describing 
microstates by the numbers of quanta {uq, ni,n2, ■ ■ ■ ,ni, . . .) = n in each single-particle 
momentum state these states being discrete if the gas is imagined to be confined to a 
finite volume V, and ordered by their kinetic energies so that pf < Pj if i < j- The total 
number of particles N ~ Ui, and the number density v — N/V. We presume that 
there is a unique state \pq) of zero kinetic energy (the ground state). The condensate, if 
present, will be characterized by the occupation number uq of the ground single-particle 
state being macroscopic, i.e. comparable to N. The kinetic energy \p^\'^/2m of all other 
states is higher than that of the ground state by at least 0{H'^ /2mV^^^); while this is small 
for macroscopic volume V, it is larger than kuT/N by a factor that is 0(iV^/'^) at fixed 
V = N/V. Therefore, a macroscopic occupation number in any state other than |pq) is 
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always exponentially suppressed. 

The energy of a general microstate n is 



i>0 i^j 



y- 



i>0 



2m 



y- 



i>0 



2m 



V 
2V 



[2N^ ~ nl] 



(5) 



The double sum in the first line represents the exchange energy for pairs in distinct states, 
and is rewritten as n^)^ — in the second line. Now 
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Following the remarks in the last paragraph, nj is microscopic for all j > 0. So the error 
in the last line of (|5]) grows more slowly with N than E{n, 1/) as TV — >^ oo at constant 
N/V and therefore can be neglected in the thermodynamic limit. This step is essential, 
because the negative sign in front of nf would otherwise prevent us from summing 
independently over all the with a chemical potential fi to enforce N — N, as one 
does to evaluate the grand-canonical partition function of a noninteracting gas, because 
the sums would diverge. But with the approximation (|5]), one can eliminate N and no 
in favor of the intensive variables t'o(M, t) = fio/V and v{^,, T) = N/V, leaving V as 
the sole extensive variable, and thereby obtain the grand-canonical partition function of 
the interacting system, and hence the pressure. Here i^o < is the number density of 
condensate particles only, not to be confused with the number density of the background 
state in the wave-drag calculation of the previous section. 




Figure 1. Scaled pressure P versus scaled number density v for interaction strength 
e = 1 [eq. j6}]. 



One finds that the condensate appears when the number density is greater than i^crit (T) = 
C{'i/2){mk^T/hhfl'^, which is the same as thecritical density as for a noninteracting 
gas. In terms of scaled number densities v = v/vcrit and i>o = ^'o/^'criti ^ scaled pressure 
P = P/vciitksT cx P/T^^'^, and a dimensionless measure of the repulsion term 

0^-^^„it =iC(3/2)4^, (6) 
ksT ^ AdB 
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the equation of state is 

P = _ ij>2) + Li5/2(z)/C(3/2), (7a) 
i>= Li3/2(z)/C(3/2) +max(i>o,0), (7b) 
z — exp(— 6'i>o) if and only if i>o > 0. (7c) 

Here Lis(z) is the polylogarithm [Lis(z) = J2kLo / ^'^ where this converges], which 
relates to the Bose-Einstein integrals when s is an integer or half-integer The parameter 
X is a sort of fugacity. When ~ 0, meaning that the condensate is absent, then z 
is an implicit function of i) defined by eq. ( fTa] ); otherwise, eq. ( |7c| i gives z explicitly. 
The noninteracting gas is recovered by setting 6 = 0: then pressure is constant once the 
condensate appears, if temperature is fixed. 

A plot of the equation of state (|7]i is shown in Fig.[T]for 9 = 1, corresponding to strong 
scattering, although the steps leading to eqs. (|7]i probably cannot be justified unless 6 1. 
The condensate is absent along the red section, and present along the blue. The part of the 
curve beneath the dashed line segment is unphysical, and the endpoints of this segment 
define two phases in contact; it can be shown that <f PAv = around the loop defined 
by the dashed segment and the unphysical lobe below it. The density jump at the phase 
transition scales cx 6 when 6* ^ 1. In that limit, however, then because of the factors 
of in eqs. ( |7b] l and (jTcJl, the pressure rises only slowly with increasing v > i^crit until 
1/ > 9~^, at which point // w i^g because z and the polylogarithms become small. 

In other words, while the actual phase transition is rather weak for weak collisionality 
(Cscatt ^ -^Ib)' ™ '■h^ same limit an isothermal system goes from non-degenerate to 
almost entirely degenerate over a small increase in pressure but a large increase in density. 
It is this feature of the equation of state that leads to constraints on the collisionality of 
RDM from galactic rotation curves. 

6. RDM halos 

The simplest model for dark halo in virial equilibrium that resembles inferences from ob- 
servations of galaxies is an isothermal sphere. For a classical ideal gas or system of colli- 
sionless WIMPs, the isothermal sphere has a unique density profile, apart from rescalings 
of the central density and core radius. It is the limit of the sequence of Emden spheres, 
as is also discussed in Chandrasekhar's first book ||3T|, and is characterized by constant 
velocity dispersion and constant P/ p at all radii. Although the density profile cannot be 
obtained in closed form, p{r) ^ r^^ at r 3> this implies a mass contained with radius 
r that is asymptotically linear, M(r) ^ r, and therefore a rotation velocity for circular 
orbits, Vc{r) = GM{r)/r, that is asymptotically constant, in general agreement with ro- 
tation curves of spiral galaxies out to the largest observable radii. The total mass profile 
(DM-nbaryons) appears to be remarkably close to p cx r^^ within at least the luminous 
regions of elliptical galaxies ll34l . 

Using the equation of state (|7|, one can construct self-gravitating isothermal spheres of 
RDM. The condensate dominates within the core radius Tc = tto, and the density profile 
there is approximately the same as for zero temperature, namely an ti = 1 Emden model. 
At much larger radii, p ex. as for the classical isothermal sphere. The behavior at 
intermediate radii is quite different from the classical case and depends strongly upon the 
collisionality parameter (|6]). This is shown in Figure |2] All of these models are scaled 
to the same central density, core radius, and asymptotic rotation velocity, which turns out 
to require that the scaled central number density i>(0) cx 9^^; so apart from dimensional 
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Figure 2. Self-gravitating isothermal spheres of RDM for various collisionalities 6 
[eq. l|6j] . Solid curves and left axis: mass density. Dashed curves and righthand axis: 
rotation curves. 

scalings, they form a one-parameter sequence. Paradoxically, the strongly collisional 
model, 9 — 1, looks most like the classical case, but for 6* ^ 1, the density drops sharply 
at the edge of the core and is approximately constant out to r ss The explanation 

for this behavior is that the doubling of the repulsion due to exchange forces acting on the 
nondegenerate component within the core (the nondegenerate component is never entirely 
absent) renders this component almost unbound, so that its density becomes small and 
declines slowly until such a radius that its own interior mass begins to be larger than that 
of the core. This structure is analogous to that of a red giant, in which the central parts 
are supported by degeneracy pressure and there is a large increase in entropy across the 
hydrogen burning shell, leading to a distended convective envelope. 

The drop in density at the edge of the core is associated with a drop in the rotation 
curve, which eventually recovers and is approximately constant beyond O^^I^Tc- Fig. [2] 
shows that drops by a factor ^ 2 for 6 = 10^"'. A drop this large or larger would seem 
to be in contradiction with what is inferred for the contribution of the halo to after 
that of the baryonic matter is subtracted. Thus, insofar as RDM halos are (or would be) 
isothermal, we can safely conclude that 6 > 10^'^. 

7. Limits on the RDM boson mass 

Eliminating Uq between (Tscattt — {mUo)"^ /t^It-^ and Vc = ira = ttUq / AGm^ leads to 
a = (4/7r^)G^fi.~*m^r(!. Therefore, one can derive an upper bound on m in terms of the 
constraint inferred by ||28l for ct/to from the Bullet Cluster: 

m<7xl0-^( "^"l '^'eV/c'. (8a) 

On the other hand, taking AdB ~ h/mvc in eq. (|6|, where Vc is the asymptotic circular 
velocity of the halo, one has a lower bound on m in terms of the constraint 6 > 10^'^ 
obtained in the last section: 

This corresponds to AdB S 0-05 cm for Vc = 100 km s^^. 




10 



Chance and Chandra 



The two inequalities ( |8a] l & ( (8b] l are clearly incompatible for any astronomically ac- 
ceptable values of core radius and circular velocity. Therefore, it seems that repulsive 
dark matter, at least of the sort we have been considering here, can be ruled out. 

The lower bound ( [8b] l derives from our assumption that the halo is isothermal, however, 
so one should ask whether that assumption is reasonable. Local thermodynamic equilib- 
rium should be established on a collision time icoii ^ Kufp/vc, but global isothermality 
is achieved on a conduction time tcond ^ r"^ / {vc\-nip)- If the mean-free path were small 
compared to galactic scales, one could have tcond ^ ^coii- As noted in our discussion 
earlier, however, cr/m w 1 gcm^^ corresponds to A^fp r in the solar neighborhood, 
so if eq. ( [Sa] ) were an approximate equality, then these timescales would be comparable 
to one another, and rather less than a Hubble time. On the other hand, if the mass were 
very much less than ( [Sa] ), then the RDM would suffer little scattering, so that even tcoii 
could be longer than the age of the Galaxy. Thermodynamics considerations would not 
apply. And yet, despite the lack of scattering, the repulsion could still provide a mini- 
mal Tc ^ 1 kpc since (as we have already noted) the core radius and the scattering cross 
section involve different combinations of Uq and m. 
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